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摘 要 : 风沙 流 的 风 程 效应 研究 是 定量 获取 风沙 流 沿 程 变化 的 核心 和 难点 , 风 程 效 应 是 指 输 沙 率 随 


或 地 块 长 度 的 增加 而 不 断 增 大 ,而 后 趋 于 稳定 的 变化 特征 ,饱和 输 沙 率 (f,) 和 饱和 小径 


长 度 (有 是 风 程 效应 的 重要 参数 。 采 用 自动 连续 称 重 式 集 沙 仪 , 以 河北 圾 上 地 区 康 保 县 境内 典型 
旱 作 农田 为 研究 对 象 , 观 测 了 2017、2018 年 和 2021 年 内 4 次 典型 风蚀 事件 ,分 析 近 地 表 5S cm 高 度 
风沙 流 的 风 程 效应 在 Smin 时 间 尺 度 下 的 变化 特征 。 结 果 表 明 :(1) 近 地 表 输 沙 通 量 随 风 程 距离 的 
增 大 而 增 大 。(2) 4 次 风蚀 事件 中 到 ,的 变化 范围 在 11~280m 之 间 , 并 存在 明显 差异 ,其 变化 与 风速 


无 关 。(3) 近 地 表 风沙 流 的 扣 . 与 风速 (UV) 呈 老 函 数 关系 。(4) 风 程 效应 的 变化 特征 与 地 表 可 蚀 性 因 


子 、 地 表 微 地 魏 变 化 有 着 紧密 联系 ,未 来 应 对 不 同 的 土壤 类 型 和 质地 农田 的 风 程 效应 进行 深入 
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风蚀 是 造成 土壤 退化 和 粉尘 排放 的 重要 原因 
之 一 ,准确 理解 风沙 流 的 变化 规律 是 风蚀 定量 研 
究 的 重要 环节 。 风 沙 流 从 可 蚀 边界 开 始 , 沿 着 下 风 
向 运 移 , 输 沙 通 量 随 着 距离 的 变化 而 变化 的 过 程 为 
风 程 效应 ” , 输 沙 通 量 随 着 下 风向 运输 距离 的 增 大 
而 增 大 ,直到 输 沙 率 达 到 一 个 稳定 值 ,也 就 是 饱和 
输 沙 率 (f,,)。 输 沙 率 从 0 到 最 大 平衡 饱和 态 之 间 ， 
风 所 运动 的 距离 为 饱和 路 径 长 度 (L,)”。 早 在 20 
世纪 30 年 代 ,Chepil 等 所 就 对 风 程 效应 进行 描述 : 风 
的 输 沙 能 力 随 着 运输 距离 的 增 大 而 增 大 。 人 们 通 
过 时 外 试验 、 风 洞 试验 和 数值 模拟 等 方法 ,对 不 同 
环境 下 风 程 效应 的 机 理 .模式 和 影响 因素 进行 了 大 
量 研 究 , 已 取得 了 很 大 的 进展 。Chepil 提 出 “雪崩 
机 制 ?解释 了 几 十 米 高 空 的 风 程 效应 。 在 这 个 理论 
中 , 跃 移 的 风蚀 颗粒 撞击 地 表 , 击 溅 起 更 多 的 运 移 


收 稿 日 期 : 2021-12-13; ”修订 日 期 : 2022-03-21 


颗粒 。 在 风速 稳定 的 情况 下 , 输 沙 通 量 随 风 程 距离 
的 增加 而 呈 指 数 增长 ,直到 风沙 流 达到 饱和 状态 。 
Owen “提出 了 一 种 空气 动力 学 反馈 机 制 , 基 于 大 量 
风 洞 试验 验证 , 风 在 对 地 表 作用 过 程 中 ,将 其 动能 
转移 到 跃 移 颗粒 中 ,增加 了 地 表 的 空气 动力 学 粗粮 
度 。 由 于 较 高 的 空气 动力 学 粗糙 度 ,导致 更 大 的 动 
能 传递 回 地表 ,进一步 增加 了 风沙 流 的 运输 能 

Paltz 等 "引入 物理 学 中 “ 弛 瑰 过 程 "的 概念 来 解释 
这 种 机 制 。Gillette 等 ”在 Owen 研究 的 基础 上 ,通过 
田间 试验 ,结合 实时 气象 数据 .得 沙 通 量 以 及 表层 
土壤 性 质 , 提 出 了 土壤 阻力 机 制 来 定量 地 解释 更 大 
空间 尺度 上 的 风 程 效 应 。 风 程 效 应 模型 应 用 于 定 
量 模拟 研究 不 同 机 制 下 的 风沙 流 。 初 期 风 程 效应 
的 探讨 通常 基于 假定 的 理想 状态 下 ,其 中 的 参数 为 
风速 .气象 条 件 结 皮 、 团 聚 体 粒 径 等 ”“"。 许 多 研 
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究 通过 风 洞 实验 “"” .数值 模拟 等 方法 验证 ,这 些 模 
型 在 描述 自然 环境 中 风 程 效应 的 机 制 和 影响 因素 
时 存在 较 大 误差 "”。Bauer 等 基于 海岸 海滩 构建 
了 三 角 葡 数 模 型 ,模拟 狭窄 海滩 风沙 流 的 风 程 效 
应 ,但 仅 适 用 于 海滩 宽度 小 于 六 ,的 狭 罕 海 滩 环 境 。 
Stout 玉 根据 风沙 流 发 育 的 自 平衡 原理 提出 了 指数 
生长 模型 以 描述 某 一 高 度 上 颗粒 物 沿 沙 床 表面 运 
移 过 程 , 能 够 较 好 地 模拟 风沙 流 沿 程 变 化 的 “ 快 - 
慢 " 发 展 模式 。 但 其 忽视 了 风沙 流 中 颗粒 碰撞 以 及 
颗粒 撞击 地 表 对 风 程 效应 产生 的 影响 ,因此 难以 描 
述 自 然 环 境 下 的 沿 程 风沙 流 "”。Andreotti" 在 Stout 
模型 的 基础 上 进行 了 修正 ,并 通过 风 洞 模拟 和 田间 
测量 2 种 方法 来 研究 LZ, 与 风速 的 关系 ,结果 表明 ， 
LL 的 变化 和 风速 无 关 。Fryrear 等 "提出 指数 增长 
模型 ,用 于 描述 不 同 下 风 距 的 风沙 流 , 输 沙 率 沿 输 
沙 距离 的 增 大 呈 “ 慢 - 快 - 慢 ” 的 发 展 模 式 。 

近 地 表 风沙 流 运 输 机 制 受 多 方面 因素 的 影 
响 。Lynch 等 ”在 北 爱尔兰 海滩 进行 研究 ,对 风 的 
负 有 反馈 机 制 进行 描述 。 实 验 表明 输 沙 通 量 在 短 距 
离 内 增长 迅速 , 随 着 地 块 长 度 逐 渐 增 加 , 输 沙 通 量 
稳定 衰减 。Zhang 等 在 腾 格 里 沙漠 东南 部 对 L, 进 
行 研究 ,发现 输 沙 通 量 随 着 地 块 长 度 的 增加 , 先 增 
大 .后 减 小 再 增 大 。Zobeck 等 和 Stout 等 ”结合 
地 表 因 子 和 气象 数据 进行 实地 研究 ,发 现 可 蚀 性 土 
壤 干 团聚 体 颗 粒 和 地 表 空 气动 力学 粗糙 度 、 以 及 测 
量 高 度 是 L, 的 主要 影响 因素 。 在 风 洞 试验 中 ,有 充 
足 沙 源 供给 , 干 松 土壤 特性 2 、 持 续 吹 风 时 间 和 风 
速 大 小 是 L, 的 主导 因素 。 在 盐湖 和 戈壁 ,地 表 湿 
度 、 土 壤 物 理 ( 生 物 ) 结 皮 也 会 影响 L,”*”。 在 海岸 
沙滩 地 带 ,海滩 几何 形状 和 风速 对 LL 起 着 主导 作 
用 下。 这 些 机 制 揭 示 了 地 表 可 人 蚀 性 因子 状况 与 风 
因子 状况 决定 着 风 程 效应 的 发 展 变化 过 程 。 


近 地 表 风沙 流 的 精准 定量 观测 仍 是 一 个 迫切 
需要 解决 的 问题 。 以 往 的 研究 中 ,无 论 是 经 过 多 次 
改进 的 Bagnold 集 沙 仪 ”, 还 是 国际 广泛 应 用 的 
MWAC 集 沙 仪 中 和 BSNE 集 沙 仪 * 都 属于 机 械 型 的 
集 沙 仪 ,必须 在 一 场 风 蚀 事件 或 一 段 时 间 后 ,人 工 
收集 风蚀 物 样品 , 带 回 实验 室 称 量 才能 得 出 风沙 流 
观测 数据 下。 受 仪 右 水 平 限 制 , 风 程 效 应 的 研究 通 
稼 是 基于 较 大 的 时 间 尺 度 ,一 些 风 程 效 应 主要 研究 
的 相关 时 间 尺 度 如 表 1 所 示 : 

研究 发 现 ,采用 较 大 的 时 间 尺 度 研 究 风 程 效 
应 ,一些 关键 变化 特征 难以 避免 地 因 平 均 化 处 理 而 
被 忽略 ? ,野外 环境 中 风沙 流 的 平衡 饱和 态 是 不 断 
变化 的 , 风 程 效应 研究 受 时 间 尺 度 的 影响 比较 显 
H ,选取 较 小 的 时 间 尺 度 进行 观测 能 够 更 加 详细 地 
捕捉 风 程 效应 参数 变化 特征 。 基 于 此 ,本 研究 使 用 
高 分 辩 率 自动 连续 称 重 式 集 沙 仪 呈 ,分 析 近 地 表 5 
cm 高 度 的 风沙 流 $ min 时 间 尺 度 的 时 空 变化 特征 ， 
并 对 风 程 效应 的 影响 因素 进行 讨论 。 


1 材料 与 方法 


1.1 研究 区 概况 

河北 坝 上 地 区 属于 北方 农 牧 交错 带 ,地势 相对 
较 高 ,全 境 平 均 海拔 1450 m, 位 于 中 温带 亚 干旱 区 ， 
属 典 型 温带 大 陆 性 气候 。 年 均 温 1.2 % ,年 均 降水 
量 350 mm, FRE 1772 mm ,无 霜 期 114 d, 年 均 风 
3E 2.99 ms ', 是 风蚀 的 易 发 区 ,主要 发 生 在 春季 |。 
该 区 域 土地 利用 方式 以 耕地 林地 .草地 为 主 ,是 传 
统 农 业 区 与 畜牧 业 区 耦合 作用 区 ,农田 风蚀 灾害 
严重 ,是 京 津 刺 地 区 重要 沙 源 地 之 一 。 农 田地 表土 
壤 层 蒲 , 质 地 松散 且 较 为 均一 ,以 沙 质 土 和 沙壤土 
为 主 ,风蚀 强烈 。 本 研究 于 2017.2018 年 和 2021 


表 1 风 程 效应 主要 研究 的 相关 时 间 尺 度 


Tab. 1 Main studied related temporal scale for monitoring the fetch effect of aeolian sand transport 


试验 地 点 集 沙 仪 类 型 时 间 尺 度 参考 文献 
得 克 萨 斯 州 ,美国 BSNE 单 次 风蚀 事件 Fryrear 等 
fe SEB TIN , E BSNE 单 次 风蚀 事件 Gregory 等 ™ 
得 克 萨 斯 州 ,美国 BSNE 单 次 风蚀 事件 Hagen ^ 
格林 威 治 沙丘 ,爱德华 王子 岛国 家 森林 公园 ,加拿大 Vertical tzopsi? T Bauer 4^ 
腾 格 里 沙漠 ,中 华人 民 共 和 国 EAI 10 min, 20 min Zhang 45 
河北 省 康 保 县 ,中 华人 民 共 和 国 平 口 集 沙 仪 单 次 风蚀 事件 牛 艳 频 ” 
萨 顷 托 ,阿根廷 MWAC 单 次 风蚀 事件 Sterk Ai 
马 吉 里 根 海岸 , 北 爱 尔 兰 电子 称 重 集 沙 仪 “ 15 min Lynch 等" 
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年 观测 了 4 次 风蚀 事件 ,观测 场地 位 于 河北 坝 上 地 
区 康 保 县 境内 典型 旱 作 农田 ,场地 1 位 于 六 段 村 
(41°55'37"N , 114°27'28"E, ,海拔 1480 m) ,场地 2 位 
于 马场 (42*07'54"N ,114°48'22"E, HK 1303 m) (A 
1a) ,具体 观测 时 间 及 地 点 如 表 2 所 示 。 在 后 文中 ， 
将 2017 年 4 月 19 日 试验 称 为 风蚀 事件 1,2018 年 4 
月 30 日 试验 称 为 风蚀 事件 2,2021 年 4 月 15 日 试验 
称 为 风蚀 事件 3,2021 年 5 月 4 日 试验 称 为 风蚀 事件 4。 

按照 当地 耕作 习惯 ,试验 场地 为 旱地 ,于 每 年 3 
月 下 旬 进 行 翻 耕 把 平 ,实验 时 段 内 地 表 均 无 植被 覆 
被 ,是 典型 翻 耕 农田 ,4 次 试验 前 地 表 的 可 蚀 性 颗粒 


(a) 研究 区 观测 点 位 置 


图 例 
A 观测 场地 
土地 利用 类 型 


(b) 事件 1 仪器 布设 


图 例 


( <0.84 mm) 分 别 占 80.26% 、78.28% 、83.71% 、 
70.19% ,说 明 4 次 风蚀 事件 的 地 表 均 具有 高 度 可 蚀 
性 。 试 验 场地 大 致 为 西北 一 东南 走向 。6 段 试验 场 
地 的 上 风向 为 天 然 草地 ,北面 为 小 麦田 ,下 风向 为 
耕地 。 马 场 试验 场地 的 西北 Ri A EAE F 
风向 为 耕地 。 为 了 避免 灌 丛 沙 堆 影 响 迎 风气 流 ,在 
试验 中 采用 防 尘 网 制作 的 防 尘 屏 障 建造 了 一 个 20 
m 宽 的 人 为 不 可 蚀 边 界 ,布置 在 集 沙 仪 上 风 辐 10 
m lb, 
12 实验 仪器 与 方法 

4 次 风蚀 事件 所 采用 的 集 沙 仪 均 为 河北 师范 大 


马场 


>z 


(c) 事件 2 仪器 布设 


图 例 
梯度 风速 仪 


(g) CCST 


注 :CCST(Cyclone type continuous-weigh sand trap ) 为 自动 连续 称 重 式 集 沙 仪 ;Met one 为 三 维 超声 风速 仪 。 
图 1 试验 样 地 位 置 及 仪器 布设 


Fig. 1 Location of the test sites and layout of the test instruments 
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表 2 风蚀 事件 基本 信息 
Tab.2 Basic information of wind erosion events 


风蚀 事件 时 间 地 点 


事件 1 2017 年 4 月 19 H 10:40—14:10 六 段 村 (场地 1) 
事件 2 2018 年 4 月 30 日 12:10 一 15:45 马场 (场地 2) 
事件 3 2021 年 4 月 15 日 13:25 一 14:55 马场 (场地 2) 
事件 4 2021 年 5 月 4 日 14:45 一 16:15 马场 (场地 2) 


学 自主 设计 人 研发 的 自动 连续 称 重 式 集 沙 仪 (Cy- 
clone type continuous-weigh sand trap, CCST) , 如 图 
le.g 所 示 , 集 沙 仪 由 连续 称 重 模 块 .数据 采集 模块 
和 电池 模块 构成 。 称 重 模 块 最 大 量程 为 400 g, 分 辨 
率 0.01 g,CR6 数 据 采集 器 (Campbell Scientific Inc.) 
采集 频率 可 调节 ,最 高 达 50 Hz ,野外 试验 采用 1 Hz 
的 采集 频率 ,能够 实现 对 风沙 流 的 高 频 连 续 观 测 。 
CCST 集 沙 口 朝向 可 以 根据 风向 自行 调整 ,试验 中 6 
个 集 沙 口 均 沿 主 风 向 设置 ,高 度 均 设置 为 0.05 m. 
风蚀 事件 1 和 2 中 的 CCST 集 沙 仪 均 布 设 在 距 不 可 
蚀 边 界 下 风向 20 m、40 m、70 m、100 m、130 m 和 
160m 处 ,风蚀 事件 1 的 集 沙 口 朝向 设 为 270" ,风蚀 
事件 2 的 集 沙 口 朝 向 设 为 310"。 风 蚀 事件 3 和 4 的 
CCST 集 沙 仪 布设 在 距 不 可 蚀 边 界 下 风向 10-110 m 
处 ,每 个 集 沙 单元 间距 20 m, 以 4x6 和 矩阵 布设 , 集 沙 
口 朝向 设 为 315"。 由 于 CCST 集 沙 仪 在 使 用 中 可 能 
HARV OHE .被 风 歇 牌 或 数 采 失 灵 情 况 ,因此 ， 
在 使 用 中 将 选取 数据 质量 较 好 的 一 列 进行 分 析 。 

野外 观测 风速 采用 2 种 风速 仪 ,风蚀 事件 1 和 2 
的 观测 采用 北京 师范 大 学 研制 的 便携式 梯度 风速 
仪 …, 如 图 所 示 ,9 个 风 杯 高 度 设置 为 0.05 m ,0.1 
m,0.3 m,0.5 m,0.8 m,1.0 m,2.0 m,3.0 m #11 4.0 m, 
风速 分 辨 率 为 0.1 ms ,并 采用 CR1000 数 据 采集 器 
记录 1 min 平 均 风速 数据 。 梯 度 风 速 仪 计算 方法 
如 下 : 
将 实测 风速 数据 通过 最 小 二 乘法 进行 拟 合 ,得 
到 计算 风速 的 公式 : 

U,=A+BlnZ (1) 

式 中 :凡是 高 度 为 Z 时 的 风速 (ms 1);4、B 为 回归 

风蚀 事件 3 和 4 的 观测 采用 三 维 超声 风速 仪 
(Met one) ^" ,如 图 1h 所 示 , 测 量 高 度 为 1.2 m, 测 量 
频率 为 50 Hz, 数 据 采 集 频 率 为 1s。 仪 右 均 布设 于 
距 不 可 蚀 边 界 60 m 的 位 置 。Met one 风速 仪 计算 风 
速 的 过 程 如 下 : 


由 于 电子 仪 需 在 野外 应 用 中 存在 不 稳定 现象 ， 

在 数据 分 析 前 ,首先 计算 出 1 min 的 平均 值 ,并 将 无 
效 值 吻 除 ,然后 根据 水 平 风 速 计算 出 风向 偏 角 : 

0=arctan (u,/v,) (2) 


a 7 arctan (u,/w.) (3) 
式 中 ;页 .页 、 而 为 实测 坐标 系 xwy.z 下 的 平均 三 维 
风速 时 间 序 列 wm's” );6 为 垂直 偶 角 ;ao 为 水 平 侦 角 。 
将 Met one 风 速 仪 坐 标 根据 偏 角 调 整 至 来 流 反 
向 与 主 风向 一 致 ,进而 得 到 新 的 坐标 系 , 计 算 过 程 
如 下 : 


ui=u,cosO+v, sind (4) 
v, 7 -u, sinÜ +v, cos 0 (5) 
w =w, (6) 
u =u cosa +v sing (7) 
v=v (8) 


式 中 :ww 为 xy 坐标 轴 调 整 偏转 角度 后 的 三 维 
KE (m-s) ;wv,、w: 为 实测 坐标 系 下 的 三 维 风 速 
(mss) ;uw 为 x.z 坐 标 轴 调 整 偏转 角度 后 的 三 维 风 
Hi (ms), 

以 1 min 为 时 间 间 隔 ,计算 主 风 加 风速 值 CU， 
m:s ) ,计算 公式 如 下 : 

U= NE +7 (9) 

IP: g, 五 为 坐标 轴 调 整 偏转 角度 后 水 平和 垂直 方 
向 风速 的 平均 值 (ms"')。 

Fryrear 和 Saleh 提 出 的 风 程 效 应 模型 适用 于 模 
拟 实际 环境 下 的 风沙 流 沿 程 变化 ,模型 表达 了 输 沙 
通 量 的 沿 程 变化 遵循 慢 - 快 - 慢 的 变化 规律 , 呈 S 型 
曲线 增长 ,并 且 其 关键 路 径 长 度 不 同 "。 风 程 效 应 
模拟 方程 如 下 : 


FO) =F! -e*"| (10) 
SOE fa) Py EES x ERRUR (g * em? min!) ;fi 为 
饱和 输 沙 率 (g: em? min) , fra 0.63/(x) ; b NRE 
地 块 长 度 (m)。 

应 用 RR 语言 软件 的 最 小 二 乘 非 线 性 回归 
(Gauss-Newton 迭代 算法 ) 计 算 5 min 时 间 尺 度 下 4 
次 风蚀 事件 的 输 沙 通 量 , 得 到 f/,, 和 6, 并 采用 决定 系 
数 ( 尼 ) 和 均 方 根 误差 (RMSE) 检 验 模型 的 拟 合 效 
果 。 饱 和 路 径 长 度 (人 是 风 程 效应 模拟 中 的 一 个 
重要 参数 , 指 当 输 沙 通 量 接近 饱和 时 的 地 块 长 度 ， 
即 当 所 -f(x)<f/100 时 相应 的 地 块 长 度 "。 
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2 结果 与 分 析 


2.1 近 地 表 风沙 流风 程 效应 特征 

分 析 发 现 , 输 沙 通 量 随 着 集 沙 时 间 的 推移 ,有 
逐渐 增 大 的 趋势 (图 2)。 风 蚀 事件 1.2.4 在 初始 1~ 
2 左右 , 输 沙 通 量 增 长 速度 达到 最 大 ,然后 逐渐 放 
慢 , 能 够 较 好 地 反映 风沙 流 输 沙 率 沿 程 的 变化 规 
律 。 输 沙 通 量 增幅 大 的 时 刻 与 风速 增 大 的 时 刻 基 
本 对 应 ,说明 风速 是 影响 风沙 流 搬运 能 力 的 重要 因 
素 。 输 沙 通 量 随 风 程 距离 的 增加 而 增加 ,但 风蚀 事 
件 270 m 处 累积 输 沙 通 量 小 于 100 m 人 处 ,这 可 能 是 由 
于 100 m 处 集 沙 口 朝向 在 大 风 条 件 下 发 生 仿 移 , 影 
响 集 沙 效 率 。 风 蚀 事件 3 试验 场地 11:30 开 始 下 雪 ， 
地 表 仿 水量 增 大 ,根据 风速 仪 和 集 沙 仪 数据 反映 ， 


— 风速 
(a) 风蚀 事件 1( 风 速 ) . - ”平均 风速 13.36 mis? 


E 
E: 
区 
6 1 1 1 1 1 
10:30 11:40 12:50 14:00 15:10 
时 刻 (北京 时 间 ) 
(b) 风蚀 事件 1( 累 积 输 沙 通 量 ) 
+ 200 
8 一 20m 
& 150+ — 40m 
hu: — 70m 
E) 100 | — 100m 
S — 130m 
= 50 — 160m 
E 
Bk 0 
10:30 11:40 12:50 14:00 15:10 
时 刻 (北京 时 间 ) 
- 20 . ( 风蚀 事件 3( 风 速 ) — 风速 
--- 平均 风速 10.89 ms” 
E 
区 10 
11:40 12:50 14:00 15:10 
时 刻 (北京 时 间 ) 
< do ( 风蚀 事件 3( 累 积 输 沙 通 量 ) 
B — 30m 
2 — 50m 
at — 70m 
A 50 + — 90m p 
= — 110m 
= ea 
BR 0 I I 1 
11:40 12:50 14:00 15:10 
时 刻 (北京 时 间 ) 


在 12:00 一 14:00 虽然 风速 较 大 ,但 输 沙 通 量 几 乎 没 
有 增长 ;14:00 后 , 随 着 风速 的 增 大 , 累积 输 沙 通 量 波 
动 增长 。 

分 析 2017 年 4 月 19 日 10:40 一 14:10,2018 年 4 月 30 
日 12:10 一 15:45,2021 年 4 月 15 H 13:25—14:55, 2021 
年 5 月 4 日 14:45 一 16:15 的 沿 程 输 沙 率 数据 (图 3~5) 
发 现 ,5 min 时 间 尺 度 内 , 近 地 表 输 沙 率 大 体 上 随 着 
地 块 长 度 的 增加 而 增加 。 并 且 , 输 沙 率 增长 呈 波 动 
上 升 的 趋势 ,总 体 上 沿 程 变化 特征 为 " 慢 - 快 - 慢 ? 曲 
线 ,大 致 呈现 出 S$ 型 曲线 的 增长 规律 ,这 与 Fryrear 等 ” 
提出 的 风 程 效应 模型 相 吻 合 。 
2.2 饱和 输 沙 率 与 风速 的 关系 

通过 风 程 效应 模型 的 预测 能 够 较 精确 地 反映 
该 农田 地 块 的 风沙 流 沿 程 变化 特征 。 利 用 R 语言 软 


风蚀 事件 2( 风 速 ) 一 风速 
(c) 风蚀 事件 2( 风 速 ) 一 平均 风速 11.22 ms 


E 
a 12 
区 6 
13:00 14:00 15:00 16:00 
时 刻 (北京 时 间 ) 
z g (d) 风蚀 事件 2( 累 积 输 沙 通 量 ) 
9 — 20m — 130m 
> — 40m — 160 
= 60 二 70m 
图 40 一 100m 
z 
= 20 
BK 0 i meer a 1 j 
13:00 14:00 15:00 16:00 
时 刻 (北京 时 间 ) 
«39 。 (加 风 亿 事件 4( 风 速 ) 一 风速 
--- 平均 风速 11.88 ms 
14:10 15:00 15:50 16:40 
时 刻 (北京 时 间 ) 


(b) 风蚀 事件 4( 累 积 输 沙 通 量 ) 
一 10m — 90m 


E 

b 20 - — 30m — 110m 

a — 50m 

tl 15H 

10 

m s 

BR 0 í 1 1 1 
14:10 15:00 15:50 16:40 

时 刻 (北京 时 间 ) 


图 2 风速 和 不 同 距 离 处 累积 输 沙 通 量 


Fig.2 Wind velocity and cumulative sediment flux at different distances 
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图 3 风蚀 事件 1 输 沙 率 的 沿 程 变化 


Fig. 3 Variations of aeolian transport rate along the path of wind erosion event 1 
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图 4 风蚀 事件 2 输 沙 率 的 沿 程 变化 


Fig.4 Variations of aeolian transport rate along the path of wind erosion event 2 
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(d) 风蚀 事件 4(14:45 一 15:15) 
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注 : 因 风 蚀 事 从 
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F 3 中 距 不 可 蚀 边 界 10 m Ib RUA AS TP DL fe ous] 


和 件 中 堵塞 ,所 以 将 10 m 人 处 数据 删除 。 


图 5 风 人 名 事件 3 和 风蚀 事件 4 输 沙 率 沿 程 变 化 


Fig. 5 Variations of aeolian transport rate along the paths of wind erosion event 3 and event 4 


件 中 的 非 线性 回归 分 析 对 4 次 风蚀 事件 5 min 时 间 ”其 他 3 次 的 模拟 效果 一 般 。 人 研究 发 现 , 在 风蚀 事件 


尺度 的 近 地 表 .进行 预测 ,采用 Fryrear 模 型 "进行 
BWE, J HEHA EI RA RMSE 检验 模型 的 拟 


REIFE „fan Lea BADER , EXT foros AYE 


合 情 况 。 经 过 相关 性 检验 后 ,将 RR<0.40 B5 Zi 


BRO AP, 


EES 和 事件 4 由 于 仪器 原因 ,只 取 第 3 d 


列 数据 进行 计算 。 计 算 结果 如 表 3 所 示 : 


RR 越 接近 1, 拟 合 效 果 越 好 ,RMSE 越 接 近 0, 误 
ZEB), AR R^«0.40 的 数据 后 ,Fryrear 模 型 预测 结 


RHEI RH 0.86, FI RMSE 为 0.25。 模 型 对 于 ” 现 , 近 地 
风蚀 事件 2 的 风 程 效应 有 着 较 好 的 模拟 效果 ,对 于 风蚀 物 运 移 主 要 分 为 2 个 阶段 ,初始 阶段 为 增长 阶 


程 变化 有 显著 的 影响 。 


论 


3.1. 近 地 表 风沙 流风 程 效应 的 时 间 尺 度 
通过 分 析 4 次 风蚀 事件 中 风沙 流 的 治 程 变化 发 


表 累 积 得 沙 通 量 随 着 风速 的 增 大 而 增 大 ， 
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X3 5 min 时 间 尺 度 的 模型 拟 合 结果 及 统计 学 参数 特征 


Tab.3 Model fitting results and statistical parameter characteristics at 5 min time scale 


风蚀 事件 fag cm” min” b/m Lalm R RMSE 
事件 1 最 大 值 5.97 68.099 104.00 0.84 1.35 
平均 值 3.01 51.099 78.93 0.70 0.71 
最 小 值 0.57 43.604 67.00 0.58 0.16 
事件 2 最 大 值 1.18 169.800 258.00 0.96 0.72 
平均 值 0.56 57.480 87.61 0.81 0.18 
最 小 值 0.12 26.510 41.00 0.54 0.01 
事件 3 最 大 值 18.26 184.300 280.00 0.97 1.10 
平均 值 5.17 91.870 90.18 0.85 0.35 
最 小 值 0.07 18.970 58.51 0.67 0.01 
事件 4 最 大 值 2.02 95.360 127.00 0.86 0.15 
平均 值 0.34 27.250 39.54 0.62 0.05 
最 小 值 0.08 7.090 11.00 0.40 0.01 


段 ,以 流体 侵蚀 为 主 , 随 着 地 块 长 度 的 增加 , 输 沙 率 
迅速 增 大 ;第 二 阶段 为 弛 驳 阶 段 ,流体 侵蚀 率 随 着 
下 风 距 离 的 增 大 而 减 小 , 跃 移 沙 粒 撞击 侵蚀 率 逐 渐 
增 大 到 极 大 值 , 输 沙 率 达 到 饱和 状态 ,地表 风蚀 物 
运输 逐渐 达到 平衡 ,风蚀 作用 显著 ,风沙 流 在 这 一 阶 
BOIS fous, Loar FE: UD Th DR EUER B gc BP BL”, 
Jf- E. ,根据 集 沙 仪 数据 显示 , 随 着 地 块 长 度 的 增加 ， 
输 沙 通 量 有 逐渐 增 大 的 趋势 ,这 与 Taylor AY OF 
究 结果 相 一 致 。 

在 研究 中 发 现 ,5 min 时 间 尺 度 下 , 近 地 表 的 输 
沙 率 沿 程 变化 呈现 出 明显 的 波动 性 , 风 程 效应 显 
著 。4 次 风蚀 事件 的 输 沙 率 均 表现 出 随 着 地 块 长 度 
增加 而 增长 ,在 达到 峰值 后 出 现 降低 趋势 的 特征 。 
风蚀 事件 1 的 输 沙 率 沿 程 变化 呈现 出 “ 慢 - 快 - 慢 ? 的 
增长 趋势 (图 3) ,近似 凸 型 。30 m 距 离 内 , 输 沙 率 增 
长 较 慢 ,30~130 mm 处 迅速 增长 达到 峰值 ,并且 风沙 


注 : 训 为 饱和 输 沙 率 ;为 伯 和 路 径 长 度 记 为 关键 地 块 长 度 ; 尼 为 决定 系数 ;RMSE 为 均 方 根 误差 。 


水 平 ,根据 模型 的 计算 结果 表明 , ALBIS TR] VJ f... HS 
变化 范围 在 5.21~18.26 g- cm? * min” Zz [B] , La AEE 
范围 在 96~280 m Z [8] ,风沙 流 搬运 能 力 非 党 强烈 ， 
风蚀 作用 显著 。 风 蚀 事件 4 的 输 沙 率 呈 现 出 * 快 - 
慢 - 快 ” 的 增长 趋势 (图 5$d~f) ,符合 S$ 型 增长 曲线 ,在 
70m 范 围 内 , 输 沙 率 保 持 稳 定 增长 ,70~90 m 的 输 沙 
率 有 所 降低 ,90~110 m 的 输 沙 率 增长 速度 较 快 , 输 
沙 率 始终 保持 较 低 水 平 。 人 研究 发 现 ,5 min 时 间 尺 度 
下 ,fa 明显 波动 。 风 速 是 影响 输 沙 率 的 重要 因素 ， 
拟 合 回 归 分 析 4 次 风蚀 事件 中 的 所 .及 其 影响 因素 
风速 之 间 的 关系 发 现 ,f 随 着 风速 的 增加 而 增加 
(图 6)。 根 据 统计 学 中 展 的 计算 结果 显示 ,fi 与 风 
速 之 间 呈 现 出 显著 的 过 函数 关系 ,并 且 时 间 尺 度 越 
大 ,其 相关 性 越 强 。 

根据 模型 计算 结果 表明 , Fryrear BI X9 Br pe 
型 增长 趋势 的 风沙 流 有 着 较 好 地 模拟 效果 。 但 在 


流 达到 饱和 状态 ,130~160 mm 距离 内 的 输 沙 率 明 显 
降低 ,风沙 流 重新 趋 于 稳定 ;10:40 一 12:10 输 沙 率 表 
现 为 直线 增长 ,12:10 一 14:10 输 沙 率 增长 过 程 表 现 
出 轻微 的 波动 。 风 蚀 事 件 2 的 输 沙 率 表现 出 “ 快 - 
慢 - 快 - 慢 ”的 阶梯 型 增长 趋势 (图 4) ,波动 性 更 加 显 
著 ; 输 沙 率 在 70 m 处 达到 第 一 次 小 峰值 ,随后 降低 ， 
在 100~130 mm 处 输 沙 率 迅 速 增 长 ,直到 第 2 次 更 高 
的 峰值 出 现 ,130~160 m 处 输 沙 率 降 低 。 风 蚀 事件 3 
在 13:25 一 14:25 的 变化 过 程 与 第 一 次 相似 (图 
5a-c) ,14:25 一 14:55 输 沙 率 在 测量 的 120 m 范围 内 
呈现 出 迅速 增长 的 趋势 ,并 且 输 沙 率 始终 保持 较 高 


风蚀 事件 4 中 , 输 沙 率 在 80 mm 内 基本 达到 饱和 状 
态 ,实测 得 沙 率 仍 呈现 持续 增长 的 趋势 。 并 且 根据 
表 1 ,风蚀 事 件 4 的 尼 为 0.62,RMSE 值 为 0.05。 统 计 
"E EYF, RRE 1, RMSE 越 接近 0, 预测 值 与 实测 
值 的 拟 合 效果 越 好 ,这 表明 Fryrear 模 型 可 能 无 法 
很 好 地 模拟 5 min 时 间 尺 度 下 S 型 增长 趋势 的 风沙 
流 。 风 沙 流 沿 程 变化 过 程 受 风速 .地表 粗糙 度 和 士 
壤 性 质 影响 。5 min 时 间 尺 度 观测 研究 发 现 , 风 沙 流 
的 风 程 效应 明显 ,ZL 波动 性 较 强 ,除了 受 风 速 的 影 
响 , 地 表土 壤 干 团聚 体 构 成 .土壤 含水 率 、 地 表 粗 烽 
度 等 因素 的 变化 对 风 程 效应 的 影响 也 不 可 忽视 。 
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图 6 饱和 输 沙 率 与 风速 的 震 函 数 拟 合 关系 


Fig. 6 Power function relationship between the saturated aeolian transport rate and wind velocity 


然而 ,目前 对 风蚀 事件 过 程 中 地 表 风 蚀 影 响 因子 的 
观测 方法 尚 不 成 熟 , 尤 其 是 土壤 干 团聚 体 构成 和 地 
表 粗 烟 度 ( 微 地 貌 ) 的 观测 尚 处 于 探索 阶段 站。 
此 ,未 来 农田 地 类 风 程 效应 研究 中 ,应 着 重 关 注 地 
表 风 蚀 影 响 因子 的 实时 变化 对 风 程 效应 的 影响 。 
3.2 Lt 的 影响 因子 

对 工 , 影 响 因 素 的 探究 ,多 数 是 基于 有 定常 风 
XE .地 表 平 坦 , 颗 粒 物 松散 的 可 控 风 洞 环境 下 得 出 
的 ,L 的 变化 机 理 在 自然 环境 下 进行 探究 一 直 是 一 
个 吸 待 解决 的 问题 。 自 然 环 境 中 土壤 表面 特征 存 
在 多 样 性 ,以 及 风 的 汗 流 状况 ,导致 风 程 效 应 的 发 
展 变化 通常 表现 得 非常 复杂 。 出 于 对 土壤 表面 状 
况 实时 监测 难度 较 大 的 限制 中 ,目前 进行 实时 定量 
的 L, 的 研究 存在 不 足 。 本 研究 基于 农田 地 类 的 4 
次 风蚀 事件 中 ,5 min 时 间 尺 度 下 近 地 表 风沙 流 的 
有 ,进行 探究 。 数 据 显 示 ,上 .有 着 显著 的 变化 趋势 ， 
波动 起 伏 较 大 ,如 表 1 所 示 ,变化 范围 在 11~280 m 


之 间 。 将 其 与 可 能 的 影响 因素 风速 进行 探讨 ,发 现 
其 发 展 变化 不 依赖 于 风速 。 

在 风 因 子 状 况 不 稳定 的 野外 环境 下 , Stout 577 
认为 ,地 表 特 性 和 风力 因子 的 变化 影响 着 1h 时 间 
尺度 下 风 程 效应 的 变化 规律 。Lynch p] 15 min 
时 间 尺 度 对 海岸 海滩 进行 研究 ,研究 表明 海岸 海滩 
的 L, 主 要 受 沙 滩 表 面 湿 度 、 风 速 和 海滩 几何 形状 的 
影响 。 笔 者 认为 ,在 农田 条 件 下 , 近 地 表 的 L, 受 多 
方面 因素 的 影响 ,包括 风 因子 、 跃 移 颗 粒 位 移 距离 
和 地 表 状 况 (土壤 干 团聚 体 粒 径 组 成 .地 表 湿 度 、 田 
间 芍 高 和 地 表 植 被 履 被 等 ) 等 ,但 最 终 起 决定 性 因 
素 的 是 其 中 最 慢 的 限制 性 因子 5。 在 风蚀 发 生 过 
程 中 ,土壤 表面 性 状 也 在 实时 发 生变 化 ,直接 影响 
着 大 气 清流 ,L, 也 因此 而 发 生 改 变 。 在 本 人 研究 中 ， 
由 于 仪器 及 观测 技术 的 限制 ,对 地 表 性 状 的 监测 尚 
存在 不 足 , 因 此 ,在 后 续 的 研究 中 ,要 对 地 表 可 蚀 性 
因子 进行 更 多 的 研究 ,在 风蚀 事件 发 生 过 程 中 , 实 
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时 地 测量 农田 土壤 表层 的 变化 。 此 外 ,在 对 风 程 效 
应 的 研究 中 发 现 , 仪 絮 精度 对 实验 有 着 很 大 的 影 
响 。 本 次 实验 采用 的 CCST 集 沙 仪 属于 电子 类 集 沙 
仪 ,优点 是 数据 采集 频 度 高 ,能 够 实现 实时 连续 地 
监测 输 沙 通 量 的 变化 ,但 在 大 风 条 件 下 会 发 生 拌 
动 ,产生 震动 误差 ” ,在 数据 计算 前 需 进 行 平 均 化 
处 理 ,六 是 输 沙 通 量 达 到 饱和 时 的 路 径 长 度 ,一 些 
较 大 的 输 沙 通 量 数据 可 能 会 在 数据 处 理 中 被 过 滤 , 
从 而 影响 和 Ls 的 计算 结果 。 总 结 以 上 ,在 研究 
中 要 重视 时 间 变 异性 ,对 近 地 表 的 农田 尺度 上 对 各 
类 影响 因子 展开 讨论 ,建立 风沙 流 的 三 维 模型 ,对 
多 种 的 土壤 类 型 进行 更 丰富 的 研究 ,以 揭示 风蚀 农 
田地 表 六 ,的 影响 因素 与 沿 程 变化 规律 。 


4 结 论 

本 研究 基于 河北 坝 上 农田 地 类 2017—2021 年 
的 4 次 典型 风蚀 事件 ,对 近 地 表 5 cm 高 度 风 沙 流 沿 
程 变 化 进行 分 析 , 并 应 用 Fryrear 模 型 模拟 近 地 表 的 
/a 和 L。 主 要 结论 如 下 : 

(1) 近 地 表 的 风沙 流 输 沙 通 量 受 下 风 癌 距离 的 
影响 显著 。 初 始 阶 段 输 沙 通 量 随 着 地 块 长 度 的 增 
加 迅速 增长 , 当 风 程 距离 达到 一 定 长 度 时 ,风沙 流 
达到 稳定 状态 , 输 沙 通 量 逐 渐 达 到 饱和 ,增长 速率 
BNE, 风蚀 作用 显著 。 

(2) 4 次 事件 风 程 效应 有 着 明显 差异 , 治 程 输 沙 
率 呈 现 出 凸 型 .阶梯 型 及 S 型 3 种 主要 变化 趋势 。 
根据 统计 学 相关 系数 尺 和 RMSE 值 表明 ,Fryrear 模 
型 能 够 较 好 地 拟 合 输 沙 率 呈现 阶梯 型 和 同型 变化 
趋势 的 风 程 效应 ,对 于 S$ 型 曲线 的 风 程 效应 的 拟 合 
效果 较 差 。 

(3) 风蚀 事件 1、2、3、4 中 的 变化 范围 分 别 为 
0.57~5.97 g cm ^ * min ', 0.12-1.18 g* cm ^* min ', 
0.07-18.26 g- cm ?^* min ,0.08-2.02 g- cm?-* min ,分 
析 4 次 风蚀 事件 的 观测 结果 ,揭示 了 5 min BERI EE 
FS us ERI DER ,f 与 风速 呈 窒 函数 关系 ,并 
且 时 间 尺 度 越 大 ,两 者 间 的 究 函 数 相 关 关 系 越 强 。 

(4) 风蚀 事件 1.2、3、4 中 工 ,的 变化 范围 分 别 为 
67~104 m、41~258 m、58.51~280 m,11-127 m。4 次 
风蚀 事件 中 工 , 的 发 展 变化 具有 明显 波动 性 ,并 且 均 
不 依赖 于 风速 ,说 明 在 野外 环境 中 ,L 的 实时 定量 
变化 是 多 种 因子 共同 影响 的 结果 ,其 变化 与 地 表 可 


刨 性 因子 及 地 表 微 地 貌 的 的 实时 变化 有 着 紧密 联 
Ro 

(5) 在 野外 非 定 常 风 场 环境 中 ,LL 的 发 生机 制 
更 为 复杂 ,其 影响 因素 有 待 进 一 步 研究 。 综 上 所 
YE ,需要 在 丰富 的 环境 中 对 不 同 土壤 类 型 和 质地 进 
行 更 深入 的 研究 ,进一步 探究 L, 的 影响 因素 及 其 发 
展 变化 机 理 。 
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Abstract: The fetch effect in aeolian sand transport is a cornerstone for wind erosion modeling and sand dune 
evaluation. The fetch effect refers to an increase in aeolian sand transport with downstream distance over an 
eroding surface. The maximum sand flux (fmx) and corresponding critical fetch length (Co are the primary 
parameters related to the fetch effect. In this paper, a cyclone type instantaneous weighing aeolian sand trap 
(CCST) was used to observe four typical wind erosion events in 2017, 2018, and 2021 in two farmland regions of 
Kangbao County, Hebei Province, China. We analyzed the 5 min characteristics of the fetch effect of sand flux at 
a height of 5 cm above the ground surface. The results revealed that: (1) The aeolian sand flux near the surface 
increases with the increase in fetch length. (2) The variations of L. for the four wind erosion events range from 
11 m to 280 m and were obvious and independent on wind velocity (U). (3) The relationship between fix and U 
followed a power-law relationship. (4) The extent of the fetch effect may be related to factors including the soil 
wind erodibility and surface microrelief. More studies are required to explore how the factors related to soil wind 
erodibility and surface microrelief affect the fetch effect. 


Key words: soil wind erosion; fetch effect; maximum sand flux; critical fetch length 


